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The hadron loop effect is conjectured to be important in understanding discrepancies between 
the observed states in experiments and the theoretical expectations of non-relativistic potential 
model. We present that, in an easily operable procedure, the hadron loop effect could shift the 
poles downwards to reduce the differences and provide better descriptions of both the masses and 
the total widths, at least, of the radial quantum number n = 1 charmed and charmed-strange 
states. The l^Pi - 1^ Pi mixing phenomena could be naturally explained due to their couplings 
with common channels. The newly observed D states are also addressed, but there are still some 
problems remaining unclear. 
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I. INTRODUCTION 



Discoveries of more and more charmed or charmed-strange states in experiments attract great interest on the 
theoretical side, because some members of them have unexpected properties. In the Particle Data Group (PDG) 
tabled), six lower charmed states, £1°, D*(2007)°, D^{2400)°, Di(2420)°, Di(2430)°, D|(2460)°, and their partners 
have already been established. Recently, some evidences of three new charmed states, -D(2550), -D(2610), and Z)(2760) 
have been reported by the BABAR Collaboration H, whose features lead to intense discussions and theoretical 
suggestions of the further experimental investigations [3|-|8| . There are also nine charmed-strange states quoted in the 
PDG table among which some states' quantum numbers are undetermined. The mass spectra of these charmed and 
charmed-strange states are roughly depicted in the predictions of the non-relativistic potential model in the classic 
work by Godfrey and Isgur (referred as GI in the following) Q. However, the observed masses are generally lower 
than the predicted ones. For example, the biggest discrepancies happening in both spectra are the I^Pq states. The 
£'0(2318) is about 80MeV lower than the expectation, while the D*q{2317) is about 160MeV lower. There are a body 
of theoretical efforts at solving this problem usually by changing the representation of the potential (see Ref. \10h 
[mI and their references). Lattice calculation has also been made to explain the experiments [l^ [l^. However, the 
present systematic uncertainty of the Lattice calculations does not allow determinations of the charmed mesons with 
a precision less than several hundred MeV. 

Another expectation to shed light on this problem is to take the coupled channel effects (or called hadron loop 
effects) into account, which plays an important role in understanding the enigmatic light scalar spectrum and their 
decays [13, . In the light scalar spectrum, the strong attraction of opened or closed channels may dramatically shift 
the poles of the bare states to different Riemann sheets attached to the physical region and the poles on unphysical 
Riemann sheets appear as peaks or just humps of the modulus of scattering amplitudes in the experimental data. 
The mass shifts induced by the intermediate hadron loops have also been testified to present a better description of 
the charmonium states |19l - [23 ]. The coupled channel effects have already brought some insights into the nature of 
the charmed-strange Dsj (2317) and some other states [23ri26j . However, although this effect could explain some of 
the observed charmed or charmed-strange states, there is still a concern that this effect may also exist in those states 
previously consistent with the theoretical expectation [27| . In this paper, we will address this point by considering 
the mass shifts, induced by hadron loops, of all the firmly established charmed and charmed-strange states. Here 
we propose an easily operable way, in which we use the imaginary part of the self-energy funtion calculated from 
the quark pair creation (QPG) model }28H30l | in the dispersion relation to obtain the analytically continued inverse 
propagator and extract the physical mass and width parameters, and then apply it to the charmed and charmed- 
strange spectra to interpret their masses and total decay widths in a consistent way. It is found that the results of 
their masses and total widths are consistent with the experimental values, at least for the non-radially-excited states. 
This picture gives a natural explanation to the l^Pi — I'^Pi mixing by coupling with the same channels instead of 
using a phenomenological mixing angle. This scheme has some similarities to the methods used by Heikkila etal. [Toj 
and Pennington et al. in their study of the charmonium and bottomonium states, but there are significant differences 
with them, as discussed in the text. 

The paper is organized as follows: In Section II, the main scheme and how to model the decay channels are briefly 
introduced. The mixing mechanism is introduced in Section HI. Numerical procedures and results are discussed in 
Section IV. Section V is devoted to our conclusions and further discussions. 

II. THE SCHEME 

We start by considering a simple model at the hadron level, in which the inverse meson propagator, P^^(s), could 
be represented as [l^ 



P'i(s) = ml - s + H(s) =ml-s + J2 ^n{s), (1) 

n 

where uiq is the mass of the bare qq state and H„(s) is the self-energy functions for the n-th decay channels. Here, 
the sum is over all the opened channels or including nearby unopened channels ("just virtual"). H„(s) is an analytic 
function with only a right-hand cut starting from the n-th threshold Sth,n, and so one can write its real part and 
imaginary part through a dispersion relation 

ReH„(.) . r dzl^, (2) 

where V J means the principal value integration. The pole of P(s) on the unphysical Riemann sheet attached to the 
physical region specifies its mass and total width of a meson by its position on the complex s plane, usually defined 
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as Spoie = (Mp - iTp/2)'^. 

One could recover a generalization of the familiar Brcit-Wigner representation, usually used in experimental anal- 
yses, from Eq.([T]), as 

F-\s)^m{sf -s + imBw'^totis), (3) 

where m{s)'^ = rn^ + Ren(s) is the "runing squared mass" and Ttot{s) — lTaIi{s) / msw ■ rnsw is determined at the 
real axis where ra^s)"^ — s = is fulfilled. The mass and width parameters in these two definitions give similar results 
when one encounters a narrow resonance, but they differ when the resonance is board or when there are several poles 
interacting with each other. 

Based on the Cutkosky rule, the imaginary part of the self-energy function is expressed through the couplings 
between the bare state and the coupled channels. The relation could be pictorially expressed as FigHJ 
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FIG. 1. The imaginary part of the self-energy function. J dll means the integration over the phase space. 



Thus, one key ingredient of this scheme is to model the coupling vertices in the calculation of the imaginary part of 
the self-energy function. The QPC model, also known as the ^Pq model in the literature, turns out to be applicable 
in explaining the Okubo-Zweig-Iizuka (OZI) allowed strong decays of a hadron into two other hadrons, which are 
expected to be the dominant decay modes of a meson if they are allowed. It is not only because this model has proved 
to be successful but also because it could provide analytical expressions for the vertex functions, which are convenient 
for extracting the shifted poles in our scheme. Furthermore, the vertex functions have exponential factors which give 
a natural cutoff to the dispersion relation and we need not to choose one by hand as in Ref. [20j . 

Here, we just make a brief review of the main results of the QPC model used in our calculation. (For a more 
complete review, see [slHs^) In the QPC model, the meson (with a quark qi and an anti-quark 52) decay occurs by 
producing a quark (93) and anti-quark ((74) pair from the vaccum. In the non-relativistic limit, the transition operator 
can be represented as 

T - -37^(lml - m|00) / dVsci V4<J' (P3 + pl)3^r(^^^)xf-™0^'wo''4(P3)4(pl), (4) 

m 

where 7 is a dimensionless model parameter and 3^™(p) = p^Y[^{9p, 4>p) is a solid harmonic that gives the momentum- 
space distribution of the created pair. Here the spins and relative orbital angular momentum of the created quark 
and anti-quark (referred to by subscripts 3 and 4, respectively) are combined to give the pair the overall J^'-^ — 0++ 
quantum numbers. (/)q* — {uu + dd + ss)/^/3 and — Sij, where i and j are the SU(3)-color indices of the created 
quark and anti-quark. Xi-m ^ triplet of spin. 

Define the S matrix for the meson decay A BC as 

{BC\S\A) ^ I -2m5{Ef - Ei){BC\T\A), (5) 

and then 

{BC\T\A) = 5^{Pf - P,)m'^'j^^^^-'bMjc . (6) 

The amplitude turns out to be 

_^M,^M,^M,^(p) = jVSEaEbEc {LaMl^SaMsJJaMjJ 

Mlj^ ,Msj^ .Ml^ ,Msb '"I 

x{LBML^SBMs^\JBMj^){LcMLcScMsc\JcMj,:){lml - m|00) 
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The spatial integral I^'l^ '^iL (^) given by 

B ' C J B _)_ c /i2 + /-i3 

where we have taken P = Pb = —Pc and fii is the mass of the j-th quark. iPuaLaMlj^ (^a) is the relative wave 
function of the quarks in meson A in the momentum space. 

The recoupling of the spin matrix element can be written, in terms of the Wigner's 9-j symbol, as (sij 

( 1/2 1/2 Sc ] 

X J2iScMsaSBMsJSMs}{SMs\SAMsA;l,~rn}^ 1/2 1/2 5b \. (9) 

S,Ms [ Sa 1 S j 

The flavor matrix element is 



IJ3 



h 


h 


Ic 


h 


h 


Ib 


Ia 





Ia 



(10) 



where l2. Is, I4) is the isospin of the quark qi. 

The imaginary part of the self-energy function in the dispersion relation, Eq.(l2]), could be expressed as 



/mn^-.Bc(.)^- „/, E \M''^'.^'''B^'^'cis)\\ (11) 

where |P(s)| is the three momentum of B and C in their center of mass frame. So, 

\P{s)\ ^ V(s - (mg + mcY){s - (tob - toc)^) 

~ 2s ■ ^ ^ 

Care must be taken when Eg. dTT]) is continued to the complex s plane. Since what is used in this model is only the tree 
level amplitude, there is no right hand cut for M.-'^^^'a^^^^b^^'^^c (s). Thus, the analytical continuation of the amplitude 
obeys A^(s + ie)* — M.{s — it) — M.{s + ie). The physical amplitude with loop contributions should have right hand 
cuts, and, in principle, the analytical continuation turns to be A^(s + ie)* — M{s — ie) = A^"(s + ie) by meeting the 
need of real analyticity. M{s + ie) means the amplitude on the physical Riemann sheet (the first sheet, in language 
of the analytic S matrix theory), and A^"(s + ie) means the amplitude on the unphysical Riemann sheet (the n-th 
sheet) attached with the physical region. 

With the analytical expression of the imaginary part of the coupled channel, one will be able to extract the poles 
or the Breit-Wigner parameters from the propagators by standard procedures. In principle, all hadronic channels 
should contribute to the meson mass, as considered by Heikkila et al. in studying the charmonium states [l^. Even 
all the "virtual" channels will contribute to the real parts of n(s) and renormalize the "bare" mass. Pennington et al. 
proposed that a once-subtracted dispersion relation will suppress contributions of the faraway "virtual" channels and 
make the picture simpler [20j . Since what we consider here is only the mass shifts, we could make a once-subtracted 
dispersion relation at some suitable point s = sq. It is reasonably expected that the lowest charmed state, D", as a 
bound state, has the mass defined by the potential model, uninfluenced by the effect of the hadron loops. Its mass 
then essentially defines the mass scale and thus fixes the subtracted point. So, we set the subtracted point sq = m|,o 
or So = ('Tic + 'Tiu)^ in a practical manner. The inverse of the meson propagator turns out to be 

where rupot is the bare mass of a certain meson defined in the potential model. 



III. MIXING MECHANISM 



In this scheme, all the states with the same spin-parity have interference effects and could mix with each other. For 
example, the two J =1 states of the P-wave are usually regarded as linear combinations of ^Pi and ^Pi assignments. 
Here in considering the coupled channel effect, the mixing mechanism comes from the coupling via common channels. 
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It is also believed that the 2^ Si and l^Di states mix with each other, similar to the interpretation of the charmonium 
^-(3770) state [H. 

The inverse of the propagator with two bare states mixing with each other reads 

where M^jj{s) is the mass matrix and mbare,a represents the mass parameter of the bare a state. The off-diagonal 
terms of the self-energy function is represented by the IP/ diagram for the two mixed states. The physical states 
should be determined by the meson propagator matrix after diagonalization 

Mlagis) = a{s)-'Ml,{s)a{s), (15) 

where the mixing matrix a{s) satisfies a{s)^a{s) ~ I, i.e., ct{s) is a complex orthogonal matrix since i^{s) is 
symmetric. The a{s) matrix turns to be complex when the thresholds are open. The physical poles could be 
extracted, in an equivalent way, by finding the zero points of the determinant of the inverse propagator, that is to 
solve det{F-^{s)) = 0. 



IV. NUMERICAL ANALYSES 



The bare masses of the related mesons are chosen at the values of the GFs work [9i], as also listed in Table HIl for 
comparison. As for the dimensionless parameter, 7, and the effective /3 parameters in the QPC model to characterize 
the harmonic oscillator wave functions, we choose the same values as determined from the potential in GI's work for 
self-consistency. The constituent quark masses are Mc = 1.628GeV, Mg — 0.419GeV, and M„ — 0.22GeV. 7 — 6.9 
and the values of /3s are from Ref. [H, [s^ . The physical masses concerned in the final states are the average values 
in the PDG table. The relative wave functions between the quarks in the mock-meson states are simple harmonic 
oscillator (SHO) wave functions usually used in the QPC model calculation, which brings some uncertainties into the 
calculation, as discussed later. 

There are some further explanations for the effective /3 parameters of cu states. Godfrey and Isgur have only 
presented their results of n = 1 5 and P-wave charmed states but not provided those of the Z?-wave and radial excited 
states which is needed in our discussion of the newly observed charmed states. We can only estimate the values by 
assuming that their ratios between the /? values of the cu states are similar to the ratios in the results from the other 



research groups. For example, we find the ratios of the /3 values between different charmed states in Ref. 37 1 and 
Ref.Q are almost same. Thus, the /? values of cu states used in our calculation, except those listed in Ref. 36 1, are 
P{l^Dj) = 0.44 ± 0.02GeV, P{2^So) = 0.47 ± 0.02GeV, and ^(2^81) = 0.44 ± 0.02GeV, respectively. 

The opened or nearby unopened two-body channels taken into account in our calculations are all hsted in Table HI 
Those channels with the a meson are not considered, because a meson is not regarded as a conventional qq state 
in the potential model @. It will quickly decay into two pions and the three-body decays usually present a minor 
contribution. 



Mode 


channel 1 




0+(l^Po) 1^ 


(l^Pi) 1^ 




2+(l3p2) 0'{2^So) r 


-(2='s'i) r 


(I'i?!) 3- 




0" +0" 


Dtv 


□ 


□ 






□ 


□ 


□ 


□ 




Dri 












□ 


□ 


□ 




DsK 












□ 


□ 


□ 


1- +0- 


D'tt 






□ 


□ 


□ □ 


□ 


□ 


□ 




D'r] 












□ 


□ 


□ 




d:k 














□ 


□ 


Q- + 1- 


Dp 














□ 


□ 




Doj 














□ 


□ 


0+4-0" 


DoTT 










□ 








i+(r)4-o- 


-Di(2420)7r 












□ 


□ 


□ 


1+(S) 4-0- 


_Di(2430)7r 












□ 


□ 


□ 


2+ +0' 


D2(2460)7r 












□ 


□ 


□ 



TABLE I. The channels of the charmed states considered in this paper. 
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The masses and widths are simultaneously determined, as listed in Table HIl for the charmed states, where we 
present the pole positions as well as the Breit-Wigner parameters for comparison. Remarkable improvements of the 
shifted masses of the already established charmed mesons could be found instantly. Furthermore, the total widths 
specified by twice of the imaginary part of the pole positions are also consistent in good quality with the values in 
the PDG table. 





Expt. mass 


Expt. width 


TTlBW 


r Bw 


a/Stjo/p = M — il72 


GI mass 


0"(l^S'o) 


1867 










1880 


l-(l'Si) 


2007 ±0.16 


< 2.1 


2016 


0.02 


2016 - O.Oli 


2040 


0+(l^Po) 


2318 ± 29 


267 ± 40 


2335 


233 


2275 - 125i 


2400 


"l+(liPi)" 


2422 ± 0.6 


20 ± 1.7 


2420 


16 


2410 - 7i 


2440 


"l+(l^Pi)" 


2427 ± 40 


384ti?° 


2409 


163 


2377 - 94j 


2490 


2+(l3P2) 


2462 ± 1 


43 ±3 


2453 


26 


2452 - 12i 


2500 


0-(2iSo) 


2533 ± 14(?) 


128 ± 33(?) 


2534 


25 


2533 - 12j 


2580 


l-(23Si) 


2608 ± 5(?) 


93± 19(?) 


2525 


8 


2523 - 5i 


2640 


l-(l^Di) 


2763 ± 5(?) 


61 ±9(?) 


2730 


120 


2686 - 66j 


2820 


3-(l='D3) 


(?) 


(?) 


2735 


9 


2735 - 4i 


2830 



TABLE II. Compilation of the experimental masses and the total widths(the PDG average values [l|]) of the charmed states, 
the shifted pole positions and the mass spectrum in the GI's model 0. The experimental values of 2^So, 2^51, and l^Di are 
from Ref.f3|. Here we only list the neutral charmed states. The unit is MeV. 

The D(l^So) state is a long-lived particle in the strong interaction and there is no opened strong channel, so we 
regard it to be well described as a bound state in the potential model and choose its squared mass as the subtraction 
point of the dispersion relations. 

When we calculate the mass shift of the D{l^Si), the Z?°7r" and D^tt^ threshold are both taken into account, 
because the D~^tt~ threshold is at about 2009MeV, just 2MeV higher than the observed D*(2007)''. It is a typical 
"just virtual" channel and in principle it will contribute a significant mass shift to the bare state. If the coupling to 
the Z?+7r^ threshold is excluded, the pole mass will only be shifted to about 2031MeV using this set of parameters. 

The pole of D{1^Pq) is significantly shifted down to 2275MeV, which is 125MeV down below the potential model 
prediction. The pole width is about 250MeV, which is in accordance with the experimental value within errorbar. 

The shifted pole mass favors the BABAR and Belle results [Ulli] over the FOCUS result The D{l^Pi) and 
D{l^Pi) states stay close to each other and they both have the same quantum numbers = 1"^ and similar decay 
channels. The unmixed pole positions are at ^JW(JTp^ = 2387 — z28MeV and .ys^^ap^ = 2427 — i71MeV, respectively. 
Both of the masses and the widths of the unmixed l^Pi state have large differences from the experimental values. 
It is the effect of their couplings with common channel D*7t that significantly change their pole positions to one 
narrower and the other boarder. The poles determined by the zero points of the inverse propagator matrix are at 
^S("iip^") = 2410 — i7MeV and ^S("i3pj") = 2377 — i94MeV respectively, which characterize the two observed states 
quite well. Their related Breit-Wigner parameters agree with the experimental values better. It is interesting to 
mention that in this scheme the mixing matrix, as a function of s, are complex- valued and it is not easy to find its 
relation with the mixing angle commonly used in phenomenological analyses. Here, the approximate values of the 
mixing matrix at about 2410MeV is 

, / -0.57 -h0.26i 0.87 -H0.17i\ 

aU=(2.4iGeV)^ = Q 3^ _^ Q^^^^ Q^g^ _ Q j ' (16) 

whose imaginary parts are small compared with the real parts. If one just neglect the imaginary parts and define the 
mixing matrix as usual, one obtains 

-sin6l cos6l \ / -0.57 0.87 
cos6l sin6i ) ~ \ 0.87 0.57 

which means the mixing angle 6 ~ (29° ^ 35°). It is in agreement with the value 6 ~ 35.3° obtained by considering 
the heavy quark symmetry '41]. 

As for the other charmed states which are not quoted in the PDG table, the evidences of several new charmed 
states, D{2550), £'(2610), and 15(2760), have been recently reported by the BABAR collaboration 0|. Several groups 
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have presented their tentative interpretations of the nature of these states [3|-|7[, and we make a brief summary of 
their conclusions here. £'(2550) is assigned to the 2^5*0, but its large decay width could not be explained by the QPC 
model, the chiral quark model, and the relativistic quark model, so further experimental explorations were suggested. 
Although the potential model has predicted D{2^Si) to be located at about 2640MeV, the QPC model and the chiral 
quark model also favor L)(2610) to be a mixed state of D{2^Si) and D{l^Di) to interpret its large width. There 
are conflicting opinions about the assignment of D(2760) as the heavier mixed state of D{2^Si) and D{l^Di), or as 

In our calculation, the mass shift induced by the intermediate states also reduces the pole masses of D{2^Sq), 
down to 2533MeV, but its pole width is quite narrow compared with the experimental value, as shown in Table HIl 
However, the mass of D(2'^iS'i) is shifted too much down to about 2523MeV due to many intermediate channels opened. 
Actually, the pole is even shifted down about lOOMeV below some thresholds, and its pole width is fairly small as 
well. It seems to become a quasi-bound state due to its strong coupling with the Ditt channel. Of course, there is 
still some parameter space for the effective /3 parameter and the dimensionless coupling strength parameter 7 to be 
tuned to reduce the mass shift to fit the experiment signal, because these parameters have significant uncertainties. 
In our opinion, one possible reason why the result seems to be inaccurate is the uncertainty of the SHO function 
we used to estimate the coupling vertices, which might have a larger tail than the realistic one in the high s region, 
which will contribute to the mass shift through the dispersion relation. A more realistic wave function solved from 
the linear potential model could be more reliable to describe the meson property. However, usually, this kind of wave 
function does not have an analytical representation and it can not be easily continued into the complex s plane in our 
scheme. On the other hand, it is too early to get any firm conclusion, since these states still need further experimental 
confirmations. The pole of D{l^Di) is shifted down to 2686 — i66MeV as well, whose Breit-Wigner mass is about 
2735MeV which is closer to the mass of D(2760). Unlike the l^Pi — l^Pi case, with this set of parameters, the mixing 
mechanism due to the coupling with their common channels does not change their positions much. The unmixed pole 
of D{l^Di) could also be estimated at about 2735 but its width is narrow. 



Mode 


channel 1 




l+(l'Pi) 1^ 


-(l^Pi) 


2+(l3p2) 




DK 


□ □ 






□ 


1- +0" 


D'-K 




□ 


□ 


□ 



TABLE III. The opened and nearby closed channels of the charmed-strange states considered in this paper. 



J in ^ Lj) 


Expt. mass 


Expt. width 


y'Spoic = M-il/2 


GI mass 




1968 






1980 




2112 ±0.5 


< 1.9 


2114- Oi 


2130 


0+(l^Po) 


2317 ±0.6 


< 3.8 


2358 - Oi 


2480 


"l+(l^Pi)" 


2459 ± 0.6 


< 3.5 


2470 - Oi 


2530 


"l+(l='Pi)" 


2535 ± 0.2 


< 2.3 


2508 - li 


2570 


2+(l='P2) 


2573 ± 1 


20 ±5 


2522 - 7i 


2590 



TABLE IV. Compilation of the experimental masses and the total widths(the PDG average values [i[ of the charmed strange 
states. The unit is MeV. 

The discrepancies that happen in the charmed-strange spectrum could be well addressed qualitatively, owing to 
their coupling with the opened thresholds and the nearby unopened OZI-allowed strong thresholds in Table lHll as the 
picture proposed by van Beveren and Rupp for explaining the £>*j(2317) state [l^. Some the thresholds are opened as 
a result of the isospin breaking effects, e.g., £'*j(2317) — > -Dsvr", iD*7r", whose contributions are highly suppressed by a 
factor of about {niu — md)/ {rus — (niu + md)/2) « 1/38, where the masses are the current quark masses. The coupling 
to such thresholds will contribute tiny imaginary parts of the self-energy functions, which hardly shift the mass of the 
state and only contribute to the decay widths with an order of KeV. So we completely neglect these thresholds with 
isospin breaking effects and those OZI suppressed. When we choose the value of dimensionless strength parameter 
7 = 6.9 for the charmed-strange spectrum, the mass shifts will be a little larger. We change 7 to be around 5.5 and 
obtained the shifted masses of the charmed-strange S and P-wave states, as listed in Table ITVl One could regard 
this fine tuning procedure as a "fit" , because we only want to give a qualitative description for the charmed-strange 
spectrum. Indeed the 7 parameter in the QPC model, determined by fitting to experimental decay processes, usually 
has an uncertainty of about 30% [1, H^]. Here we only list the pole positions, because they do not differ much with 
the Breit-Wigner parameters in this case as quasi-bound states or narrow states. 
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If the isospin breaking and other weak interaction channels are unopened, £>*(2112) and £>*j(2317) are the bound 
states when the bare 1^5*0 and I'^^'i states are coupled to the "just virtual" DK threshold. They show as the poles 
on the real axis of the physical Riemann sheet. It is the coupling of the bare states to the lower isospin breaking 
DsTT^ thresholds and the other weak thresholds that shift the poles to the unphysical Riemann sheets when they are 
open in reality. When the mixing of the l^Pi and I'^Pi states is not considered, they are both the bound states at 
2478MeV and 2493MeV respectively below the D*K thresholds. The mixing owing to coupling with the common 
D*K thresholds shifts the l^Pi downwards along the real s axis, and the l^Pi moves upwards and crosses D*K 
thresholds into the complex s plane of unphysical Riemann sheet. 



V. CONCLUSIONS 



In this paper, we propose a simple procedure to extract the pole positions or determine the Breit-Wigner param- 
eters of the charmed states based on the parameters in the non-relativistic potential model, by using the analytical 
representation of the QPC model to mimic the behaviors of the imaginary part of the self-energy function of the 
meson propagator. Overall improvements could be found between the pole positions or Breit-Wigner parameters and 
the well-established charmed and charmed-strange mesons. Several charmed-strange states could be regarded as the 
quasi-bound states, due to the coupling with nearby unopened OZI-allowed thresholds. In this model, the l^Pi — l^Pi 
mixing is explained by the coupling with common channels and these resultant pole masses and widths are consistent 
with the observed values. It is worth stressing that our calculation is the first one that systematically addresses such 
a broad spectrum and the decays of the members by considering the coupled channel effects, as far as we know. This 
calculation may help to improve our understanding the charmed and charmed-strange spectra. 

There are still some differences between the shifted pole positions and the parameters of the newly observed 
states. Since at the present stage the statistics of the data is still not enough to make a firm determination, further 
experimental evidences are required for a confirmation of these mesons. 
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